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Abstract 

Numerous experimental indications suggest that the Hume-Rothery mechanism plays an important role in stabilizing 
quasicrystalline and amorphous phases. However, the exponential damping of the conventional Friedel oscillations at 
the relevant, elevated temperatures T poses a severe challenge to the HR stabilization. In order to resolve this problem 
it is shown using a Feynman diagram technique that quantum correlations in the electron sea, arising from the interplay 
of Coulomb interaction and impurity scattering, can strongly enhance the Friedel oscillations in these systems even at 
elevated temperature. The resulting corrections to the Friedel potential are in agreement with available experimental 
results on amorphous HR alloys. It is proposed to include the enhancement of the Friedel amplitude derived in the 
present work into pseudopotentials through the local field factor. 

Keywords: Electronic structure; Coulomb interaction and disorder; screening; Friedel oscillations. 
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1. Introduction 



The microscopic origin of the stability of quasicrys- 
talline phases continues to be an unresolved issue. 
While at high temperatures the long-range order of 
quasicrystals and approximants may be favored over 
crystalline order because of entropic effects JlJ , strik- 
ing interrelations between the ionic and the elec- 
tronic structures [§, || especially in icosahedral (i) 
quasicrystals seem to indicate that electronic sta- 
bilization plays an important role at short and in- 
termediate length scales of of several atomic spac- 
ings: (1) The electron density of states (DOS) has 
a pronounced, structure-induced pseudogap at the 
Fermi level sp. (2) The position k p of the main 
ionic structural peak coincides with twice the Fermi 
wave number, 2kp — k p 0. (3) In dependence of 
the composition from their constituents, quasicrys- 
talline phases are only stable in small regions where 
the condition 2kp ~ k p is satisfied. (4) In the qua- 
sicrystalline phase the electrical resistivity is sub- 
stantially larger than that of each of the elemen- 
tal constituents. Similar coincidences between elec- 



tronic and ionic structures are observed in amor- 
phous noble-polyvalent alloys [|[ ||. These findings 
may be traced back to the common feature of these 
materials that the ion system exhibits concentrical, 
shell-like density correlations, where the spacing a 
between neighboring shells coincides with the Friedel 
wavelength, Xp = n/kp — 2n/k p = a. Therefore, 
it has been conjectured that the electronic Friedel 
oscillations around an arbitrary central ion give an 
important contribution to the pair potential. As a 
consequence, the total energy of the system should 
be optimized by the ions effectively being bound in 
the minima of the Friedel potential and a concommit- 
tant pseudogap formation at Ep. The importance of 
such a Hume-Rothery (HR) stabilization mechanism 
is supported by detailed theoretical studies at tem- 
perature T = both for amorphous Q and for qua- 
sicrystalline @, ||, |9[ systems. Numerical simula- 
tions [jllj show that quasicrystalline structures can 
indeed be grown by the use of pair potentials with 
appropriate repulsive (i.e. oscillatory) parts. 

However, at the relevant temperatures where qua- 
sicrystals or amorphous structures are stable (T ~ 
10 2 . . . 10 3 K), regular Friedel oscillations are expo- 



nentially damped, and the thermal ion energy is 
so large that the experimentally observed short- to 
intermediate-range concentrical ion density correla- 
tions cannot be understood on the basis of the con- 
ventional Friedel oscillations alone. 

In the present work it is shown that the interplay 
between Coulomb electron-electron interaction and 
disorder can lead to a strong enhancement as well 
as to a systematical phase shift of the Friedel oscilla- 
tions even at finite T. Both effects compare well with 
available experimental results on amorphous HR al- 
loys, and support the validity of the HR mechanism 
even at elevated temperatures. 




Fig. 1. a) Diagrammatic definition of the diffusion ver- 
tex r. b) Polarisation 11(0, g) including leading order 
quantum correction induced by disorder and interac- 
tions. Dashed lines denote electron-ion scattering T- 
matrix tt p , the wavy line with solid triangles the ef- 
fective Coulomb interaction. 



2. Model and effective interaction 



We here discuss an effective model for the electron 
motion in amorphous and quasicrystalline systems 
at T > 0. The ionic density correlations mentioned 
above constitute a scattering potential for the elec- 
trons whose scattering T-matrix t? p i s J by defini- 
tion, proportional to the static ion structure factor, 

p (x S(\k — k'\), which in turn is peaked at a mo- 
mentum transfer q = \k — k'\ = k p ~ 2kp and thus 
leads to enhanced backscattering. As has been shown 

iimno 



the latter not only generates a pseudogap 
but at the same time leads to a substantial increase 
of the electron transport or density relaxation rate, 
t -1 , over the quasiparticle decay rate, . r _1 is re- 
lated to the conductivity a by a = ne 2 T/m*. Thus, 
t -1 > r~p is a generic feature of the amorphous 
and the quasicrystalline state. In quasicrystals, when 
the conductivity is substantially reduced below the 
Drude result, we may have r _1 3> . 

In addition, at the relevant, finite temperatures 
the phase coherence of the electrons is lost on the 
length scale of the inelastic mean free path. It fol- 
lows that the electrons cannot probe the long-range 
order in a quasicrystal. Rather, they experience an 
effective potential made up of randomly placed, spa- 
tially extended scattering centers, each one charac- 
terized by the T-matrix tt p [fl2| . In such a potential 
the electronic motion is diffusive instead of ballistic, 
similar as in amorphous metals. When the electron 
coherence length is long enough (low T), the mo- 
tion may be subdiffusive with a diffusion exponent 
/3 < 1/2 H 0. However, for the present purpose 
the precise value of (3 is unimportant, and we will 
assume j3 = 1/2 (classical diffusion) in the following. 

Diffusion, as a dissipative process, is difficult to 
incorporate in an ab initio calculation. Therefore, 
we will choose a Feynman diagram technique, where 
diffusion arises by averaging over all (quasi-)random 
configurations of the system. In a diffusive elec- 
tron sea screening is inhibited, so that the effective 



Coulomb interaction, v^'{z, Z), between electrons 
with complex frequencies z and z + Z acquires a long- 
range, retarded part Ea|, 



ff _ v q T 2 (z,Z,q) 

« 1 ' >~ e RPA (Z,q) ' 



Ane 2 



(1) 



where e RPA {Z,q) = 1 + 2™ a /{ZsgnZ" + iq 2 D) is 
the disordered RPA dynamical dielectric function, 
and the diffusion vertex, defined in Fig. 1 a), is 



?{z,Z,q) = 



i/r SgnZ" 
Z+iq 2 D sgnz" 

1 



z"(z + Z)" < 
otherwise. 



(2) 



D = 1/3 VpT and " denote the diffusion constant 
and the imaginary part, respectively, and /3 = 1/2. 
The long-range nature of Vq*f is a consequence of the 
hydrodynamic (Z,q — > 0) pole of T, Eq. (2). Since 
diffusion is a classical phenomenon, guaranteed by 
particle number conservation, the form Eq. (1) of 
the effective interaction persists at finite T. 

There are two experimental indications for the dif- 
fusion model of electron transport to be valid in i- 
quasicrystals. First, note that the diffusion-enhanced 
effective Coulomb interaction Eq. (1) implies the 
well-known y/\E — ef\ dependence of the DOS in 
the pseudogap of disordered systems Q , where the 
half-integer power is characteristic for diffusion. The 
fact that a powerlaw dip in the DOS at Ef with an 
exponent very close to 1/2 has been observed in i- 
quasicrystals by tunneling measurements Jl6[ may 
be taken as an indication that the electron motion 
is indeed diffusive in these systems, and that the ef- 
fective interaction has indeed the Alt'shuler-Aronov 
form Eqs. (1), (2). Second, the diffusion model, 
based on a finite phase coherence length, explains 
why the spikiness of the DOS, predicted for ideal 
quasicrystals at T = 0, has up to now not been ob- 
served experimentally p7|. 
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3. Electron density response 

The incomplete screening of the electron-electron 
interaction in a diffusive metal may be expected to 
drastically affect the screening charge distribution 
around an ion in the electron sea as well. In or- 
der to calculate this effect on the Friedel oscilla- 
tions, we must consider the static charge density 
response x(0, q) in the vicinity of q = 2kp. It is 
given in terms of the polarization function 11(0, q) 
as x(0,q) = n(Q,g)/(l - v q IL(0,q)). The first term 
of Fig. 1 b), n(°)(0, q), corresponds to the well- 
known Lindhard function (RPA) Jig] , In this dia- 
gram r contributes only a nonsingular factor of 0(1), 
since here the effective interaction v^' enters in the 
static limit at wave numbers q ~ 2kp, where the 
diffusion vertex T is structureless. We are thus led 
to consider quantum corrections where T gives con- 
tributions with vanishing frequency and momentum 
transfer, so that the hydrodynamic transport prop- 
erties become important, although the response is 
taken at large external wave numbers q. The most 
singular contribution of this type arises from the 
quantum correction IjWft), q), shown in Fig. 1 b), 
2nd diagram. 

While diffusive density relaxation occurs in general 
for large times, t > t, in U^(0,q) it is, in addition, 
cut off for times larger than the life time r qp of the 
quasiparticles, which arc interacting via v e J^ . Thus, 
the frequency transfer in this term extends over the 
nonvanishing (see above) range tz} < \Z\ < r — 1 . It 
may be evaluated explicitly as 113, [2(J , 

n (1) (o,^-c^^x (3) 



1/(4T) 



(2Trh) 2 {s f t)V 2 
x - 1 



cosh 



2 _v_ 

IT 



l(x-iy + ( 



-)2]3/4 ' 



where x = x{v) = (q/2kp) / ' + v / 'ep and C is 
a numerical constant of 0(1). It is seen that for 
T = 0, Tgp 1 — > nW(0, q) exhibits a powerlaw diver- 
gence cx — sgn(q — 2kp)/\q — 2kp\ 1 / 2 at q = 2kp. Al- 
though at finite T or the divergence of IT- 1 ) (0, q) 
is reduced to a peak, the inverse dielectric function, 
l/e(q) = 1/(1 - v q U(0,q)), still has a q = 2k F di- 
vergence at a critical transport rate, t c _1 (T) because 
of the vanishing denominator. This leads to a sys- 
tematical enhancement as well as to a phase shift 
of the Friedel oscillations (see below). On the other 
hand, when there is no enhanced backscattering, we 
have Tqp/r — 1, and the peak structure of 11(0, q) 
vanishes. Eq. (3) constitutes an extension of previ- 
ous work |20| in that the finite quasiparticle life time 
is explicitly taken into account. The parameter r _1 
may be varied by changing the composition of the 
alloy. 
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Fig. 2. Crystallization temperature T c as a function of 
the DOS suppression at ef, (1 — g). Data points rep- 
resent T c for a,-Cui- x Sn x [5]. The solid curve is the fit 
of the present theory (see text). Vertical line: position 
of CAT. The inset shows the phase shift <p of the first 
maximum of the charge density distribution p(r). Solid 
line: theory. Data points with solid line: measurements 
[5] for a,-Cui- x Sn x - 



4. Comparison with experiments 



In the following, the theory developed in the pre- 
vious section is applied in the region l/epr qp <C 1 
to a large class of noble- polyvalent metal alloys like 
Cui-^Snj;. These HR alloys exhibit an amorphous 
to crystalline transformation (CAT) as a function of 
the composition of the alloy, and the thermal sta- 
bility may be continuously varied. Remarkably, in 
all these systems (1) the thermal stability reaches 
a maximum at or near the CAT (Fig. 2) j^], and 
(2), assuming that the ions sit in the minima of the 
Friedel potential, the measured ionic positions sug- 
gest that there is a systematical phase shift ip of the 
Friedel oscillations ||, p(r) oc cos(2kpr — ip), with 
ip = tt/2 at the CAT (inset of Fig. 2). 

Fourier transforming 1 — l/e(q) to obtain p(r) 
[Jl8| shows that for incomplete Fermi surface- Jones 
zone matching, i.e. small r" 1 ~ T ( jp 1 j t ne quantum 
corrections generate density oscillations p^ x \r) cx 
— cos(2kpr)/r 3 , which overcompensate the conven- 
tional Friedel oscillations, implying a phase shift of 
ip = 7r |2(j). As r" 1 — > t c -1 , the increasing 2k F peak 
of l/e(q) leads, in addition, to density oscillations 
p^{r) oc sin(2fci?r)/r 2 , so that in the vicinity of r" 1 



p(r) oc — 



cos(2fcpr) 
{2k F rf 



_! sin(2fc F r) 
A(t )- 



(2k F r) 



(4) 



with A^- 1 ) ~ 0.343tt(1 - r" 1 /^ 1 )- 1 / 2 . The ex- 
ponent 1/2 is characteristic for diffusive behavior. 
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Thus, the Friedel oscillations are shifted by ip = 
7r-tan" 1 [2fc F r A^ 1 )] ~ tt/2 + \/{2k F rA), i.e. the 
diverging Friedel amplitude necessarily goes hand in 
hand with tp — > tt/2 + 0. Note that, in contrast to 
the conventional Friedel oscillations, this divergence 
is robust against damping due to finite T or disorder. 
The point where the amplitude, A, diverges should 
be identified with the CAT, since at this point the 
fluctuations of the Friedel potential also diverge, al- 
lowing the system to find its crystalline ground state. 
This process explains in a natural way the observed 
composition dependence of the thermal stability and 
of the phase shift ip mentioned at the beginning of 
this section. 

For a direct comparison with experiments the con- 
trol parameter of the theory, r _1 , must be translated 
into a parameter which is experimentally accessible: 
It follows from the scattering theory jlj| ||o) that 
r _1 = t^ 1 + 7 S(2kp), where the peak of the ionic 
structure factor, S(q — 2fcp), controls the backscat- 
tcring amplitude, 7 is a constant, and r" 1 is an offset 
due to momentum independent scattering. S(2kp) in 
turn is proportional || to the measured, structure- 
induced suppression of the DOS N(sf) at the Fermi 
level, 1 — N(eF)/N (sF) = 1 — g, compared to the 
free electron gas, N (ef)- The resulting fit of the 
crystallization temperature, T c , is shown in Fig. 2, 
where the contribution to the stability coming from 
the pseudogap formation is assumed to be linear in 
(1 — g) (dashed line). The characteristic increase of 
T c at the CAT, explained by the present theory, is 
clearly seen. The inset shows the calculated phase 
shift, ip, and the measured shift of the atomic nearest 
neighbor position relative to the position of the first 
conventional Friedel minimum, a a = ir/kp. Note 
that there is no adjustable parameter in ip, once the 
fit of T c has been performed. The overall behavior 
of the shift is well explained by the theory; however, 
the experimental data approaches (p — tt/2 faster 
than predicted. The latter may be attributed to the 
fact that, as seen from the discussion after Eq. (||), 
the higher-order Friedel minima approach p = tt/2 
faster than the first one. In this light, the agreement 
between theory and experiment is remarkably good. 



5. Conclusion 

The structural similarities |Q| between amorphous 
alloys and i-quasicrystals suggest that the quantum 
effect discussed above may be important in the latter 
systems as well. In fact, quasicrystals seem to ful- 
fill all the necessary conditions for this effect to oc- 
cur, i.e. effectively diffusive electron motion |ll| [lJJ 
and t _1 3> Tqp 1 ■ The latter is supported by the 
Fermi surface matching, i.e. by the experimental 



observation [|l6) and theoretical prediction [Q, ||, || 
of structure-induced pseudogaps. Moreover, another 
more commonly known effect of disorder-enhanced 
Coulomb interaction, the \/\E — £p| dependence of 
the DOS in the pseudogap |15| , may have been al- 
ready observed in i-quasicrystals by tunneling mea- 
surements of the DOS (]]J. 

In conclusion, we have shown that the Friedel os- 
cillations are enhanced by Coulomb interaction in 
the presence of disorder and enhanced backscatter- 
ing. The results are relevant for the stability of amor- 
phous and quasicrystalline metals. It is proposed 
to include the enhanced Friedel potential calculated 
in the present work in the pseudopotential of more 
quantitative ab initio calculations. 
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P. Haussler, R. Haberkern, A. Huck, T. Kopp, and 
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supported by DFG through SP Quasikristalle. 
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